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ABSTRACT 

PID  controllers  are  used  for  decades  in  controlling  processes  in  linear  feedback  control  systems.  Their  use 
requires  accurate  and  effective  tuning  to  satisfy  an  acceptable  performance  for  the  control  system. 

Large  number  of  processes  are  classified  as  or  approximated  by  a  second  order  model.  They  may  me  under 
damped,  critically  damped  or  over  damped.  This  paper  presents  the  tuning  of  PID-controllers  used  with  second  order  over 
damped  processes.  The  process  damping  ratio  is  from  1  to  10  and  its  natural  frequency  is  from  2.5  to  15  rad/s.  The  tuning 
technique  depends  or  minimizing  the  integral  of  square  of  error  (ISE)  between  the  time  response  of  the  system  to  a  unit 
step  input  and  its  steady-state  response.  It  was  possible  to  achieve  a  PID-controlled  system  with  controllable  overshoot  and 
settling  time.  The  tuning  results  are  listed  in  tables  for  direct  use  depending  on  process  damping  ratio  and  natural 
frequency. 

The  tuning  approach  is  reduced  to  only  one  set  of  controller  parameters  to  simplify  the  tuning  process  used  with 
any  process  approximated  by  a  second-order  over  damped  model.  The  tuning  results  are  compared  with  controller  tuning 
using  standard  forms  showing  better  performance  of  the  control  system  using  the  proposed  tuning  in  the  present  work 
regarding  the  maximum  percentage  overshoot  and  settling  time  of  the  control  system  using  both  exact  and  simple  tuning 
parameters  of  the  PID-controller. 

KEYWORDS:  Over  Damped  Second  Order  Processes,  PID  Controllers,  Simple  Controller  Tuning,  MATLAB 
Optimization  Application 

INTRODUCTION 

PID-controllers  have  simple  construction  using  operational  and  new  generation  amplifiers  [1,2,3].  PID-controllers 
still  find  interest  of  control  engineering  researchers.  Sung  and  Lee  (2000)  proposed  an  identification  algorithm  for  the 
automatic  tuning  of  PID  controllers  to  guarantee  better  accuracy  and  to  provide  more  frequency  data  sets  of  the  process  [4]. 
Vrancic,  Strmcnik  and  Turicic  (2001)  used  the  multiple  integration  of  the  process  time  response  for  calculating  the 
parameters  of  the  PID  controller  [5].  Shen  (2002)  proposed  a  tuning  method  for  PID  controller  providing  the  performance 
assessment  formulas.  His  method  is  based  on  a  genetic  algorithm-based  design  technique  [6].  Tavakoli  and  Tavakoli 
(2003)  presented  a  tuning  technique  for  PID  controllers  controlling  first-order  plus  delay  processes  based  on  dimensional 
analysis  and  optimization  techniques.  They  used  ISE,  IAE  and  ITAE  performance  criteria  [7].  Gaing  (2004)  presented  a 
design  method  for  determining  the  optimal  PID  controller  parameters  of  an  AVR  system  using  the  particle  swarm 
optimization  algorithm  providing  high-quality  solution  [8]. 
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Syrcos  and  Kookes  (2005)  presented  a  general  mathematical  formulation  for  the  development  of  customized 
PID  control  tuning.  The  presented  a  number  of  case  studies  clarifying  their  proposed  methodology  [9].  Chen  and  Chang 
(2006)  studied  the  optimal  design  of  a  PID  controller  for  an  active  magnetic  bearing  using  genetic  algorithms  providing 
experimental  and  simulation  results  [10].  D'Emilia,  Marra  and  Natale  (2007)  described  a  theoretical-experimental  approach 
allowing  the  evaluation  of  the  adequateness  of  new  methods  for  auto-tuning  with  respect  to  traditional  ones  [11]. 
Ramasamy  and  Sundaramoorthy  (2008)  used  the  impulse  response  instead  of  the  step  response  of  the  plant  to  tune  the 
PID  controller,  requiring  no  approximation  of  the  plant  by  any  model. 

They  derived  formulae  for  the  calculation  of  PID  controller  tuning  parameters  [12].  Malwatker,  Sonawane  and 
Waghmane  (2009)  proposed  a  model-based  design  of  PID  controllers  for  higher-order  oscillating  systems.  They  obtained 
the  controller  parameters  from  a  reduced  third-order  model  and  presented  examples  to  demonstrate  the  effectiveness  of 
their  proposed  method  [13]. 

Pai,  Chang  and  Huang  (2010)  presented  a  simple  calculation  method  of  a  PI/PID  controller  tuning  for  integrating 
processes  with  dead-time  and  inverse  response  based  on  a  model  [14].  Matausek  and  Sekara  (2011)  proposed  a  tuning 
procedure  for  ideal  PID  controller  in  series  with  a  first-order  noise  filter  based  on  the  extension  of  Ziegler-Nichols 
frequency  domain  dynamics  of  a  process  [15].  Ayala  and  Foelh  (2012)  presented  the  design  and  tuning  of  two  PID 
controllers  using  the  non-dimensional  sorting  genetic  algorithm  approach  offering  simple  and  robust  solutions  providing 
good  reference  tracking  performance  [16].  Saglam,  Tutomo  and  Kurtulan  (2013)  presented  a  PI  controller  tuning  method 
for  cooling  of  the  hydrogen  production  unit  within  the  polymer  electrolyte  membrane  fuel  cell  based  micro-cogeneration 
system.  Their  proposed  analytical  tuning  rules  are  based  totally  on  the  model  parameters  using  the  digital  control  theory 
and  experiences  in  industrial  control  problems  [17].  Venugopal,  Ganguly  and  Singh  (2013)  used  soft  computing 
methodology  based  on  fuzzy  logic  to  design  and  tune  PID  controllers  for  better  performance  [18].  Arturo  et.  al.  (2013) 
developed  a  methodology  for  PID  controller  tuning  by  coupling  the  gain-phase  margin  method  with  the  genetic  algorithms. 

They  applied  the  proposed  technique  using  MATLAB  simulation  and  experimentally  on  CNC  machines  and  an 
industrial  robot  [19].  Baloochy  (2013)  employed  a  self-regulating  approach  using  the  relay  oscillation  method  to  specify 
process  parameters.  He  used  the  modified  Zeigler-Nichols  tuning  method  to  tune  the  PID  family  of  controllers  [20]. 
Korsane,  Yadav  and  Raut  (2014)  demonstrated  using  a  method  of  tuning  PID  controllers  using  different  tuning  techniques. 
They  implied  their  technique  to  a  PID  controller  used  with  a  first-order  plus  time  delay  [21].  Srinivas,  Lakshmi  and  Kumar 
(2014)  compared  the  performance  of  genetic  algorithms  based  PID  controller  with  conventional  PID  controller  for  various 
tuning  techniques  in  the  case  of  three  tank  level  process  for  set-point  tracking  [22].  Aborisade  and  Adewuyi  (2014)  studied 
the  temperature  control  of  a  gas-fired  oven  using  a  PID  controller.  They  implemented  the  Ziegler-  Nichols,  good  gain  and 
Skogestad's  tuning  methods  to  control  the  output  temperature  of  the  oven  [23]. 

ANALYSIS 
Process 

The  process  has  the  transfer  function: 

Mp(s)  =  cd„2  /  (s2  +  2i;<Bn  s  +  con2)  (1) 
The  second-order  process  of  Eq.l  may  be  for  an  exact  second-order  process  or  of  a  non-  oscillating  dynamic 
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system  of  order  higher  than  2  and  can  be  approximated  by  a  second-order  one.  The  damping  ratio  of  the  second  order 
system  is  in  the  range:  1  <  C,  <  10 

Controller 

The  controller  used  in  this  study  is  a  proportional  +  integral  -  derivative  (PID)-  controller. 
The  PID-controller  has  a  transfer  function  Gc(s)  given  by: 

Gc(s)  =  Kpc  +  (K/s)  +  Kds  (2) 
Where:  Kpc  =  controller  proportional  gain 
Kt  =  controller  integral 
Kd  =  controller  derivative  gain 
Control  System  Transfer  Function 

Assuming  that  the  control  system  is  a  unit  feedback  one,  its  transfer  function  becomes: 

M(s)  =  (b0s2  +  b^  +  b2)  /  (s3+  a,s2+  a2s  +  a3)  (3) 
where 

b0  =  co„2  Kd,  b !  =  co„2  Kpc 

b2  =  con2  Kt,  at=  2C,to„  +  con2  Kd 

a2  =  ran2(  1  +  Kpc),  a3  =  con2  K; 

System  Step  Response 

A  unit  step  response  is  generated  by  MATLAB  using  the  numerator  and  deniminator  of  Eq.  3  providing  the 
system  response  c(t)  as  function  of  time  [24]. 

Controller  Tuning 

The  sum  of  absolute  error  (ISE)  is  used  as  an  objective  function,  F  required  by  the  optimization  process.  Thus: 
F  =  j[c(t)-css]2(4) 

where  css  =  steady-state  response  of  the  system  (css  =  b2/a3). 

The  performance  of  the  control  system  is  judged  using  two  functional  parameters: 

The  maximum  percentage  overshoot,  OSmax 

The  settling  time,  Ts 

The  maximum  percentage  overshoot  and  the  settling  time  of  the  control  system  are  obtained  using  the 
stepinfo  command  of  MATLAB  [24].  The  optimization  procedure  used  depends  on  the  MATLAB  optimization  toolbox 
using  the  command  "fininunc".  Thus,  it  does  not  need  the  definition  of  any  functional  constraints.  In  this  case  the 
judgement  of  the  tuning  results  depends  on  the  values  of  the  step  response  maximum  overshoot  and  settling  time. 
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Tuning  Results 

A  MATLAB  code  is  written  with  process  parameters  con  and  C,  as  input  parameters  and  the  parameters  of  the 
controller  are  evaluated  by  minimizing  the  objective  function  (Eq.4)  using  the  command  "fminunc" . 
Tables  1  and  2  presents  the  results  of  the  optimal  tuning  procedure  used  in  this  work  for  process  natural  frequencies  in  the 
range:  2.5  <  con  <  15  rad/s 

Table  1:  PID-Controller  Tuning  for  over  Damped  Second-Order 
Processes  Having  a  2.5,  5,  7.5,  rad/s  Natural  Frequencies 
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Table  2:  PID-Controller  Tuning  for  over  Damped  Second-Order 
Processes  Having  a  10,12.5, 15,  rad/s  Natural  Frequencies 
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SIMPLIFYING  THE  TUNING  PARAMETERS 

Having  a  detailed  insight  to  the  controller  parameters  values  in  Tables  1  and  2  we  find  that  they  are  close  to  each 
other  except  for  the  derivative  parameter  at  some  of  the  process  natural  frequencies.  This  has  given  an  idea  of  taking  the 
average  of  all  the  controller  parameters  values.  This  step  results  in  an  average  proportional,  integral  and  derivative  gain  of: 

Kpcav  =  29.5776    with  standard  deviation  of  1.3733. 

Kmv=  31.4980     with  standard  deviation  of  1.3948.     i  (5) 

Kdav  =  31.2383     with  standard  deviation  of  1 .5576. 

COMPARISON  WITH  STANDARD  FORMS  TUNING 

To  examine  the  effectiveness  of  the  simplified  tuning  procedure  used  in  this  work,  it  has  been  compared  with  the 
results  of  controller  tuning  using  the  ITAE-based  standard  forms  [25].  The  standard  characteristic  equation  of  a 
second-order  closed-loop  control  system  of  the  form  given  in  Eq.3  is  [25]: 

s3  +  2.97co0s2  +  4.94co02  s  +  co03  =  0  (6) 
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The  PID-controller  parameters  are  tuned  using  the  standard  forms  for  a  specific  process  natural  frequency  and 
damping  ration  by  comparing  the  coefficients  of  the  characteristic  equations  in  Eqs.3  and  6,  giving: 


co0  =  (con2K1)"3 


Kpc=(4.94(B02/a)n2)-l 
Kd=  (2.97co0  -  2  i^con)  /  con2 


(7) 


Using  Tables  1  and  2  and  Eqs.5  and  7,  the  tuned  controller  parametrs  of  the  PID  controller  using  the  exact  tuning, 
simplified  tuning  and  the  standard  forms  for  a  selected  values  of  process  natural  frequency  and  damping  ratio  are  given  in 
Table  3. 


Table  3:  PID-Controller  Tuning  for  over  Damped  Second-Order  Processes 
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The  time  response  of  the  control  system  using  the  three  sets  of  PID-controller  parameters  of  Table  3  is  compared 
in  Figures  1-5: 
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Figure  1:  Unit  Step  Response  of  the  PID  Controlled  Second-Order 
Process  Having  co„  =  5  and  £  =  2 
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Figure  2:  Unit  Step  Response  of  the  PID  Controlled  Second-Order 
Process  Having  w„  =  7.5  and  £  =  1 
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Figure  3:  Unit  Step  Response  of  the  PID  Controlled  Second-Order 
Process  Having  to,,  =  10  and  £  =  2 
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Figure  4:  Unit  Step  Response  of  the  PID  Controlled  Second-Order 
Process  Having  to„  =  12.5  and  £  =  2 
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Figure  5:  Unit  Step  Response  of  the  PID  Controlled  Second-Order 
Process  Having  £2„  =  15  and  Z  =  1 

The  peformance  of  the  PID-controlled  overdamped  second-order  process  in  terms  of  the  control  system  maximum 
percentage  overshoot  and  settling  time  is  compared  in  Table  4.  The  maximum  percentage  undershoot  is  zero  for  all  the 
controller  parameters  sets. 
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Table  4:  Comprasion  of  the  Control  System  Performance 
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CONCLUSIONS 

•  It  is  possible  to  suppress  higher  oscillations  and  very  large  settling  time  associated  with  higher-order  dynamic 
processes  through  using  PID-controllers. 

•  Through  using  a  PID-controller  it  was  possible  to  control  the  time-based  specifications  of  the  closed-loop  control 
system  through  tuning  the  controller. 

•  The  integral  of  the  error  square  criterion  (ISE)  was  used  to  tune  the  PID-controller. 

•  The  tuning  procedure  was  reduced  to  the  assignement  of  a  one  set  of  controller  parameters  independent  of  the 
process  damping  ratio  and  natural  frequency. 

•  The  performance  of  the  control  system  was  compared  for  3-sets  of  tuning  paprameters:  the  exact  parmameters 
obtained  in  the  present  work,  the  simple  tuning  one  set,  and  the  tuning  parameters  obtained  using  the  standard 
forms. 

•  The  maximum  percentage  ovedrshoot  and  the  settling  time  using  the  present  tuning  parameters  were  superior  with 
respect  to  those  obtained  using  the  standard  forms  tuning  technique. 

•  The  simple  tuning  set  of  the  PID-controller  parameters  has  given  the  same  performance  parameters  as  the  exact 
sets  of  parameters. 

•  Using  the  simple  set  is  sufficient  in  tuning  PID-controllers  used  with  processes  approximated  by  an  overdamped 
second-order  model. 

•  Maximum  percentage  overshoot  was  as  low  as  0.0045  %  compared  with  about  10  %  using  the  tuning  approach  of 
standard  forms. 
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